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Contract  No.  K19628-73-C-0030 
R.F.  Gie3e 


A.  Background 

Soils  are  complex  mixtures  of  fine-grained  minerals  and  rock 
fragments  often  with  a<lmixcd  organic  matter.  The  minerals  are 
commonly  clays  with  others  which  are  resistant  to  chemical  and 
physical  weathering.  It  is  probably  true  that  the  important 
chemical  and  physical  properties  of  soils  are  due  to  the  clay  minerals 
present.  Therefore,  any  knowledge  or  understanding  of  soils 
must  ultimately  rest  on  study  of  clays  and  clay  minerals. 

Many  of  the  chemical  and  physical  properties  of  clays  and  clay 
minerals  which  set  them  apart  from  other  commonly  occurring  minerals 
result  from  a combination  of  1)  small  grain  size  (large  surface  area 
and  less  than  perfect  crystallinity),  2)  presence  of  interstitial 
water  between  mineral  grains  as  well  as  in  the  mineral  structure 
accompanied  by  various  cations,  3)  the  presence  of  hydroxyl  groups 
on  the  grain  surfaces  as  well  as  in  the  structure  and  4)  a layer* 
type  crystal  structure.  As  a result  of  variations  in  these  pro- 
perties, soils  are  variable  but  a possible  simplifying  feature  is 
that  clays  and  clay  minerals  can  be  grouped  into  a few  large 
families  each  of  which  should  have  similar  properties. 

A-l . Thermal  Blow-Off  Problem 


One  of  the  areas  where  the  unique  characteristics  of  clays  is 
of  prime  importance  is  in  the  reaction  of  soils  to  high  thermal 
fluxes  as,  for  e:  imple,  are  experienced  during  an  above  ground 
nuclear  explosion.  Approximately  one-third  of  the  weapon  output 
is  in  the  form  of  thermal  radiation  from  the  fireball.  It  is  well 
kr.ovn  that  on  many  types  of  ground  surfaces,  the  thermal  radia- 
tion results  in  the  ejection  of  dust  into  the  air.  As  a result  of 
these  radiation  induced  cuannes  in  the  s rface  and  the  first  few 
meters  above  the  ground,  the  -.hock  wave  as  it  propagates  away  from 
the  fireball  is  strongly  modified.  The  ai  narent  heating  of  the 
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B.  Objectives 

The  very  nature  of  soil  minerals  and  ospocially  the  clay 
fraction  makes  them  very  difficult  to  treat  experimentally.  An 
adequate  understanding  of  their  relative  stability  and  parti- 
cularly how  they  interact  with  water  molecules  adsorbed  on  the  sur- 
faces and  between  the  layers  is  necessary  in  order  to  interpret  and 
even  predict  experimental  results.  The  approach  adopted  during  this 
contract  has  been  that  such  an  understanding  can  be  achieved  by  a 
simple  theoretical  treatment  of  clay  minerals.  The  crystal  struc- 
ture of  most  of  the  important  types  of  clay  minerals  are  known 
although  not  always  with  the  accuracy  one  would  like.  A complete 
description  of  these  structures  requires  in  addition  a potential 
energy  function  which  allows  one  to  calculate  the  potential  energy 
given  the  atomic  positional  parameters  of  the  crystal  structure. 

A further  aim  of  the  research  performed  under  this  contract 
was  to  use  the  results  of  the  theoretical  description  of  clay 
minerals  to  explain  the  known  physical  and  chemical  properties 
of  clays  and  in  particular  their  behavior  vis-a-vis  adsorbed  water 
molecules 


B-l . Subtasks 


In  order  to  carry  out  the  task  outlined  in  the  contract  these 
subtasks  were  identified. 

Model  of  Bonding 

For  several  reasons  which  have  been  stated  in  previous  quar- 
terly reports,  it  was  decided  early  in  this  research  to  concentrate 
on  an  ionic  model  for  the  bonding.  The  greatest  justification  for 
this  is  the  simplicity  of  the  theory  and  the  ease  and  speed  with 
which  the  computations  may  be  performed  using  a digital  computer. 

The  simple  theory  is  tbat  the  atoms  are  fully  ionized  (Giese 
et  al,  Zeit.  fur  Krist.,  134,  275-284,  1971)  and  that  the  potential 
energy  of  any  such  crystalline  material  can  be  calculated  from 
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7.  ■ ‘ ionic  charge 

N - number  of  atoms  in  the  unit  coll 

A,U,D  3 constants 

f =*  unit  force  with  dimensions  kcal/inolo/angstrom 
(Busing,  Trans.  A.C.A.,  6,  57,  1970) 

The  three  terms  in  the  above  expression  are  the  coulomb 
(electrostatic)  van  der  Waals  and  repulsion  energies  respectively. 

Of  the  constants  in  those  three  terms,  only  the  Z's  are  known. 

The  others  must  be  determined  using  previously  determined  crystal 
structures.  The  process  is  very  lengthy  and  has  not  yet  been  carried 
out.  Rather,  since  the  electrostatic  energy  is  the  largest  term  the 
research  done  under  this  contract  has  involved  only  the  calculation  of 
the  electrostatic  energy. 

Computation 

These  computations  have  been  based  on  a modified  version  of  the 
computer  program  MANIOC  (Baur,  Zeit.  fur  Krist.,  Acta  Cryst.,  19, 

909,  1965) . All  subsequent  investigations  of  particular  mineral 
systems  has  used  a modular  programming  approach  i.e.  a separate 
program  has  been  written  for  each  type  of  investigation  (interlayer 
bond  energies,  polytype  stabilities,  motion  of  water  molecules,  etc.) 
which  in  effect  creates  a series  of  new  crystal  structure  representing 
intermediate  states  which  occur  during  a physical  or  chemical  change. 

The  new  structures  are  input  to  the  program  to  compute  the  potential 
energy  of  each. 

C.  Results 

Many  of  the  chemical  and  physical  properties  of  clays  with 
water,  organic  molecules  or  cations  involves  the  surface  energy 
of  the  clays.  This  is  by  definition  the  work  involved  in  separating 
the  silicate  layers  of  a clay  to  an  infinite  distance.  As  an  initial 
step  in  this  contract,  work  was  begun  on  an  algorithm  to  compute  the 
electrostatic  surface  energy  of  any  ionic  crystal.  This  involved 
calculatin'  changes  in  the  potential  energy  as  a function  of  interlayer 
distance  (for  clays  and  micas) . This  subtask  was  successfully  com- 
pleted and  the  algorithm  applied  to  various  sheet  silicate  structures 
(see  reference- •«) . 

As  an  extension  of  the  surface  energy  problem,  the  general  question 
of  the  stability  of  sheet  silicates  wan  examined  using  the  hypo  lesis 
that  a stable  crystal  structure  exists  in  some  sort  of  local  potential 
energy  minimum.  In  particular,  a stable  sheet  silicate  should  have 
1)  a net  attractive  force  between  the  silicate  layers  and  2)  have  the 
layers  superimposed  in  such  a way  that  a translation  of  one  layer  over 
the  other  in  a.y  direction  will  require  an  input  of  energy.  The 
hypothesis  has  been  tested  and  verified  for  micas,  kaolinite,  dickite, 
nacrite,  talc,  pyrophyllite  and  several  varieties  of  Al (OH) ^ (see 
refer  .nces) . 

The  role  of  hydrogen  bonding  h s in  the  past  often  been  invoked 
to  explain  the  stability  of  the  kaolin  minerals  and  yet  this  was  never 
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a very  satisfactory  hypothesis  because  the  interlayer  distances 
imply  v'ry  long  and  therefore  weak  hydrogen  bonds.  Considerable 
effort  wont  into  examining  this  question  because  this  is  central  to 
understanding  dchydroxylation  and  ultimate  dehydration  of  hydroxyl 
containing  c lays  and  micas.  The  interlayer  bonding  of  the  kaolin 
minerals  was  calculated  and  found  to  be  larger  whereas  when  F was 
substituted  for  OH  in  the  same  structures,  there  was  found  to  be 
a small  net  repulsion.  The  sample  conclusion  was  that  the  hydrogen 
ions  were  responsible  for  the  interlayer  bonding. 

Finally,  work  was  begun  on  the  study  of  a dynamic  system, 
vermiculito/water , in  which  the  adsorbed  phase  is  known  to  diffuse 
through  the  matrix.  Two  vermiculites  were  chosen,  one  with  Na  and 
the  other  with  Mg  as  the  interlayer  cations  because  1)  much  NMR 
work  has  been  done  on  the  two-layer  water  complex  and  2)  the  crystal 
structure  has  been  determined  with  reasonable  accuracy.  The  NMR 
results  showed  that  the  water  molecules  surrounding  the  Na  vermi- 
culite  rotate  with  an  activation  energy  of  8.5  kcal/mole  but  the 
nature  of  the  rotation  was  not  known.  Several  types  of  rotation 
were  considered  and  for  each  (i.e.  about  a vertical  axis,  about  a 2- 
fold  axis  of  the  water  molecule  etc.)  a number  of  models  were 
constructed.  For  example  in  treating  a rotation  about  a vertical 
axis,  the  atoms  of  the  water  molecules  were  rotated  in  10°  increments 
and  in  this  way  the  dynamic  process  of  rotation  was  modeled  by  a series 
of  static  models  each  representing  an  intermediate  step  in  the  process. 
The  potential  energy  was  calculated  for  each  and  the  change  in 
potential  energy  as  a function  of  rotation  gave  the  energy  barrier 
for  that  particular  type  of  motion.  The  one  with  the  least  activa- 
tion energy  was  for  a combined  rotation  about  a vertical  axis  (to 
(001))  along  with  a rotation  about  the  Na-0  axis.  The  calculated 
energy  is  10  kcal/mole  in  good  agreement  with  the  experimental  value. 

D.  Publications 

The  following  are  the  published  results  supported  in  whole  or 
part  by  the  contract.  (Sec  attached) 
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